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ABSTRACT 

'jj  An  adaptive  grid,  finite-volume  method  hae  been  used  to  solve  the  Navier- 
Stok.ee  equations  for  complete  (forebody  and  afterbody)  flowfielde  around  blunt 
be  lies.  The  code,  which  is  applicable  for  axisyianetric  or  two-dimensional 
flows,  allows  the  meeh  to  adjust  during  the  computation  to  provide  a closer 
spacing  of  mesh  points  in  regione  of  high  gradiente,  thus  minimizing  the  number 
of  required  computational  points.  The  solution  technique  is  explicit,  utilizing 
a maximum  time-step  advancement  at  each  grid  point  to  accelerate  convergence  to 
the  eteady  state.  The  code  has  been  fully  vectorized  for  efficient  eolution  on 
the  CYBER  203  computer.  A very  flexible  rezoning  routine  ie  ueed  to  concentrate 
mesh  pointe  anywhere  in  the  field,  either  by  a ueer-defined  weighting  function 
or  by  allowing  high  gradient  regions  to  adjust  the  grid.  The  grid  adjuetment 
routine  is  implicit  ir.  nature  and  represents  a very  email  portion  of  the  total 
computational  coet.  Currently,  the  code  rune  in  approximately  1.6  >'  f 

seconds  per  grid  point  per  iteration. 
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INTRODUCTION 

The  finite-volume  method  of  numerically  eolving  syetems  of  coneervation  laws 
has  been  successfully  applied  to  a wide  variety  of  problems  in  fluid  mechan- 
ics.^ ^ Its  ability  to  maintain  conservetion  of  mass,  momentum,  and  energy 
from  cell  to  cell,  even  in  reth.er  complex  nonorthogonal  coordinatee,  mak.ee  it 
particularly  attractive  for  use  with  adaptive  grid  techniquee.  For  the  purpoeee 
cf  this  paper,  a finite-volume  formulation  (FVF)  ie  defined  as  a diecrete 
approximation  to  a conservation  law  written  in  integral  form  which  (1)  uniquely 
defines  control  volumes  in  such  a way  that  control  volumee  (ceiie)  do  not  over- 
lap nor  are  gape  left  in  phyeical  spece  and  (2)  uniquely  definee  fluxee  and 
forces  acting  through  coll  walls  so  that  summability  without  residue7  (coneerva- 
tion) ie  guaranteed.  It  differs  from  a finite  difference  formulation  (FDF) 
only  in  th»  “.v*  a problem  is  approached.  For  example,  given  a nyetea  of  con- 
eervetion  lawe  we  might  consider  the  FDF  as  a discrete  approximation  to  the 
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differential  form  of  the  laws  while  the  FVF  is  a discrete  approximation  to 
the  integral  form  of  the  laws.  A review  of  various  FVF's  suggests  that  alJ 
FVF's  can  be  written  as  FDF's  (and  usually  appear  that  way  in  the  literature) 
but  not  all  FDF's  can  be  written  as  FVF's.  Furthermore,  discretizing  the 
divergence  form8  of  the  conservation  laws  in  general  coordinates  to  obtain  the 
FDF  does  not  guarantee  that  the  formulation  can  be  expressed  as  an  FVF  because 
conservation  is  not  maintained  due  to  some  confusion  of  how  to  properly  define 
the  metric  coefficients.  Various  ways  of  defining  the  metric  coefficients  to 
overcome  this  problem  have  been  given  in  References  5 and  9.  The  discussion 
there  leads  one  to  conclude  that  while  there  may  be  a certain  amount  of 
ambiguity  in  how  to  properly  define  or  discretize  the  metrics  the  problems  can 
be  overcome  without  a great  deal  of  difficulty.  It  should  alsc  be  noted  that 
the  FDF  of  a conservation  law  written  in  nondivergence  form  can  in  fact  be  a 
conservative  FVF.  An  example  of  such  a scheme  is  presented  in  Reference  10. 

The  discussion  up  to  this  point  has  been  concerned  with  taking  a differential 
form  of  a conservation  law,  approximating  it  by  an  FDF,  and  testing  if  it 
satisfies  the  requirements  of  an  FVF  (or  forcing  it  to  satisfy  the  requiremente 
by  suitably  defining  the  metrics) . The  approach  taken  in  the  presentation  which 
follows  is  to  start  with  the  integral  form  of  the  conservation  law,  approximate 
it  with  am  FVF  on  some  generalized  grid,  and  expand  the  formulation  to  eee  what 
FDF  results.  It  will  be  shown  that  no  special  treatment  of  metrics  is  required 
and  in  fact  no  terms  which  are  readily  recognizable  as  metric  coefficients  ever 
appear.  It  is  only  when  the  FVF  is  expanded  into  a format  more  familiar  ae  an 
FDF  do  theee  metrics  appear  as  ratios  of  the  dimensions  of  cell  wells  to  cell 
volumes.  Furthermore,  the  unexpanded  form  of  the  FVF  in  ideally  suited  for 
vectoritation. 

Within  this  framework  a grid  adaption  routine  has  been  developed  which 
greatly  facilitates  placement  of  men;  points  whera  needed.  A deecription  of 
how  the  adeption  algorithm  has  evolved  from  the  groundwork  established  in 
references  11  and  12  to  the  current  implicit  adaption  echeme  will  be  presented. 
Finally,  eample  calculations  of  supersonic  flow  over  a Viking  Aeroehell,  includ- 
ing flow  in  the  wake,  demonetrate  the  vereetility  of  both  the  FVF  end  the 
euaption  elgorithm  through  its  ability  to  concentrete  points  not  only  in  the 
boundery  layer  but  also  in  the  free  shear  layer. 

FINITE-VOIWffi  FORMULATION 

Coneider  e eystea  of  two-dimensional  coneervetion  lewe  which  can  be 
expreseed  by  the  vector  reletion 
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(la) 


(lb) 


where  Equations  (la)  and  (lb)  are  the  differential  and  integral  forms  of  the 
conservation  laws  respectively,  and  i and  j are  unit  vectore  in  the  x and 
y directions,  respectively. 

A control  volume  in  two-dimeneional  (or  axieymnetric)  space  is  determined  by 
a quadrilateral  whoee  vertices  are  defined  by  adjacent  mesh  points.  The  [i,j] 
cell  refere  to  the  cell  with  verticee  (i,j)<  (i+l#j) , (i,)^l) i (i+1, j+1) - 
Pointe  of  constant  j index  define  i coordinate  linee,  increasing  £ in 
the  direction  of  increasing  i.  Points  of  constant  i index  define  P| 
coordinate  lines,  increaeing  H in  the  direction  of  increasing  j (eee 
Fig.  1).  The  flux  through  the  wall  defined  by  the  points  <(i,j),  (i,j-*-l)> 
ie  calculated  using  . ^formation  from  the  (i,j)  vertex  on  the  predictor  step 
and  from  the  (i,j+l)  ve  tex  on  the  corrector  etep  for  odd  time  steps.  The 
order  is  revereed  for  even  time  stepe.  The  flu*  through  the  wall  defined  by 
the  pointa  <(i,j),  (i+l,j>>  follows  a similar  pattern  of  definition,  ueing 
information  from  the  (i,j)  vertex  on  the  predictor  step.  The  FVF  for  the  [i,j] 


cell  is  written: 


3U 

3t 


n*l 


i.j 


Ai.j  + F-  ^(y 


♦l,jtyi*l,j*l  " Yi+1,  j*  ' Gi*l,j(xi*l.j*l  *i*l,j' 
- Fi,j(Yi,j.l  - yi,j>  +Gi.jlxi,jn-Xi,j> 

‘ Fi,  j+l<yi+l.  j+1  ' yi.j*l>  * Gi.j>l<Xin,jU  " Xi.jn) 


+ Fi,j<Vi,j-yi,j)  " Gi.j(xin,j  - xi,j) 


(2a) 
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au 

3t 


n+1  "TTT  n+1  , „ 

. Ai,j  + Fi+l,j+l(yi+l.j+l  - yi+l,j}  ' Gi+l,j+l  Xi+l»j+l  ' 
h j 


- ‘ yi.j}  + GM+lUi^+1  ' Xi'^ 

‘ Fi+1, j+l(yi+l, j+1  ’ yi.j+l)  + Gi+1» j+l(Xi+l.j+1  i.j+l 

n+1  , . rn+l  _ x ) » o (2b) 

+ Fi+1, j yi+l, j ” yi.j  " i+1.1  i+1*3  ^ 


<]  - ■*{<.]  * ■Hi) 


(3) 


^re  Ai.  is  the  eree  of  the  [i,j]  cell.  Equetion  (2)  «u,t  epproximate  the 
integrel^orm  of  Equetion  (1).  The  firet  term  of  Equetion  (2)  Bust  epproxiraate 
the  integrel  of  3u/3t  over  the  entire  cell.  In  the  FVF  we  eet 


3u 

3t 


n+1 

i.j 


At 


ij 


(4e) 


3u 

7? 


n+1 

i.1 


n+1 


Ui+l.i+l  ~ Ui+l,j+l 


At 


i+l.j+1 


(4b) 


Note  the  uee  of  . At  which  veries  from  point  to  point  in  Equation  (4) . 

Now  coneider  e finite-difference  formulation.  PDF,  of  Equation  (1)  which 
can  be  trensformed  to 


(5) 


where  x “ x(C.H) » y 


y(C>n).  The  FDF  of  (5)  ia  written: 
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U?+J  = 0?  . - At.  . 

i»j  i.j  i<j 


(Fi+l.j  ~ r+  + ~ Gi«3)  ^ 


AC 


X.  . 

13 


AC 


+ <Fi.j+i  + <Gi,jn  - Gi,3>n  ♦ 

An  V.  Ln 


(6a) 


U"+J  - U?  . - At.  . 
1 / j * * J J 


(F.  . - F ,)n+l  <G.  . - G x ,)n+1  _ 

vj 1~l/3 r + i,j  At  1 C 

AC  Sj  ^ ySj 


.n+1  r ^n+1 

, <Fi.i  - Fi.i-i’  -,-  . '%)  - n- 

*!  n*ij  4"  .'13 


(6b) 


U11*1  - i(u"+1  + U"+1) 
l.D  2 ^ 1,3  i,3  / 


(6c) 


The  FVT  (Eqa.  (2)  and  (4))  and  the  FDF  (Eq.  (6))  are  equivalent  if 


5xij  E (yi+l,j+l  " yi+l,j’ 


<7a) 


V . - <yi, j+1  - yi+l.j*r  An/Aij 


(7b) 


Cyij  ’ <Xi+X,j  " Xi+1'j+l)  A^Aij 


(7c) 


ny..  £ Ui*l,j+1  - in/Aij 


(7d) 


(7e) 


\tj  ' <yi-1.3-l‘  yi.3-l’  iVAi-l.j-l 


<7f) 
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<Xi-l,j-l  ‘ Xi-1» j*  AC/Ai-l,j-l 


(x.  . , - x.  , . Ari/A.  , . . 

i,:-i  1-1,  :-i  1-1, j-i 


(7g> 


<7h) 


A linearized  analysis  of  Che  FDF  (Eq.  (6)),  after  it  has  been  expanded  into 
a one-step  scheme  ahowa  the  following:  (1)  Equation  (6)  is  a consiatent 
repi  ■‘sentation  of  Equation  (5),  (2)  the  truncation  error  for  variable 
At. . ia  first  order  in  time  and  can  be  expressed 


c + 


(8) 


where 


C 


3F  r 


3C 

3U  S 


3f 

3U  n» 


8g  n 
+ ny 


and 


hi;j(C,n>  - At(t,n>/At(Cij,nij> 

Note  that  E.  • 0 if  At . . - constant,  and  (3)  the  spatial  truncation  error  is 
2 A ^ J 

of  0(A  ) where  A is  proportional  to  a cell  wall  dimension  if  all  £*  , etc., 
are  f irat-order-accurate  approximations  to  etc.,  and  ♦ £^)/2,  etc.,  , 
axe  second  -order  -accorete  approximations  to  etc.  These  conditions  are 

formally  satisfied  by  Equation  (7).  For  example 


(£*  + Q/2  " Cx  + l tJyCn> £ + Jy££n!  + (9) 


where  J • x y - x y . 

£ n ^ £ 

Conaequently , it  ia  seen  that  while  the  FtJT  or  FVT  ia  second-order-accurate 
(for  At  ■ constant)  excessive  stretching  or  skewness  of  the  grid  can 
introduce  large  factors  of  second-order  terms.  Furthermore,  these  factors 
multiply  terms  like  32U/3£^,  etc.,  and  they  contribute  to  the  overall 
artificial  viscosity  of  the  method.  At  present,  grid-induced  errors  are 
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checked  by  comparing  results  obtained  for  ths  same  problsm  on  ths  different 
grids  constructed  from  mors  points  or  from  different  adjusting  parameters  (ess 
section  on  ADAPTIVE  GRID).  It  should  also  be  noted  that  while  Equation  (5)  is 
not  in  strict  conservation  form,  the  FVF  which  models  Equation  (lb)  is  con- 
servative in  ths  strict  sense  of  summation  without  residua.7  For  exanpls, 
regardless  of  the  grid  used,  the  FVF  returns  an  exact  uniform  flow  if  s uniform 
flow  is  us'd  as  an  initial  condition. 

It  is  noted  that  the  equations  and  analysis  in  this  section  era  for  two- 
diaensionel  flow.  An  axi symmetric  extension  introduces  only  minor  complica- 
tions. Finally  we  point  out  that  ths  FVF  (Eqs.  (2)  end  (4))  reduces  to 
MacCormack's  method  in  Cartesian  coordinates.  There  ere  an  infinity  of  other 
ways  to  define  calls  and  flux  through  cells  on  s general  grid,  each  of  which 
lead  to  s v lique  FVF  and  unique  definitions  of  metric  coefficients  for  that 
particular  FVF. 

BOUNDARY  CONDITIONS 

Only  a cursory  description  of  boundary  conditions  can  be  provided  due  to 
space  limitations. 

u • v ■ 0. 

Adiabatic  wall  or  constant  wall  temperature . 3p/3n  - 0. 

(n  coordinate  normal  to  wall) 

Rankine-Hugoniot  relations  for  a diacrats  moving  shock  are 
epplied.  Pressure  behind  shock  obtained  from  extrapolation 
from  interior  points. 

Limiting  form  of  differential  form  of  governing  equations 
solved  using  MacCormack's  method. 

Primitive  variables  obteined  by  axtrepolstion. 

Ths  present  grid  edjustment  scheme  has  evolved  from  ona  which  imparted 
velocities  to  grid  points  baaed  on  gradients  in  the  fiald  as  in  References  11 
and  12  to  one  which  rezones  the  computational  mesh  using  an  implicit  algorithm 
at  any  desired  frequency  (i.e.,  once  per  time  stsp  or  once  per  1000  time  stsps) . 
Originelly  svery  point  (i.j)  in  ths  computational  plans  wes  connected  to  the 
four  adjecent  points  (iil.j),  (i.jil)  by  springs  whoss  spring  constants 
Ki’  Kj  “*r®  determined  by  s function  of  tho  gredient  of  soete  dependent 
variable  between  the  points.  For  exampls,  the  spring  along  the  ith  row 
connecting  points  (i.j)  and  (i,j*l)  was  defined 


Well: 

Shock: 

Syimnetry: 
Outflow: 
ADAPTIVE  GRID 
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Kj  " 1 + clfi(j  - fi,j+ll/Ar  <10) 

where  c is  £ conetant  and  Ar  is  the  distance  between  the  points.  Grid 
points  on  the  boundariee  were  free  to  move  along  the  boundary  as  if  on  a 
frictionlese  rod.  Grid  points  on  the  comers  of  the  domain  were  held  fixed. 
Small  adjuetmente  in  the  x and  y directione  were  aseigned  to  grid  points 
as  determined  by  the  net  forcee  on  the  grid  points  in  the  reepective  directions. 
After  evsry  time  step,  the  grid  was  updated  and  checked  to  make  sure  that  no 
grid  overlap  was  imminent. 

In  a noninteractive  teet  case  where  p(x,y)  was  prescribed  to  model  the 
preseure  field  of  an  oblique  shock  crossing  a uniform  flowfield  (p(t)  = 0), 
the  grid  evolved  from  one  of  equally  epaced  rectangular  celle  to  the  one  shown 
in  Figure  2.  In  an  interactive  teet  caee  where  the  solution  p(x,y,t)  was 
allowed  to  evolve  with  time,  the  final  grid  distribution  (Fig.  3)  is  eeen  to 
be  much  more  erratic  than  that  of  the  previoue  caee.  The  preseure  field  for 
this  case  is  shown  in  Figure  4. 

By  sacrificing  the  column  to  column  (or  row  to  row)  influence  of  the  epring 
system  (i.e.,  no  springs  between  column  i and  column  i+1)  a euperior  algorithm 
in  terms  of  computational  coste  and  sraoothnees  of  grid  distribution  can  be 
constructed.  (Smoothneee  of  grid  distribution  is  a eubjective  judgment.  A 
coordinate  line  whoe?  direction  or  length  changee  erratically  from  point  to 
point  is  judged  nonsmooth.) 

Consider  the  coordinate  line  of  constant  index  i as  shown  in  Figure  5. 


(11) 


Let  f.  be  a table  of  dependent  variablee  at  the  point  e.  defined 

j ,»  J 


by 


fM  ■ -j'  fj,2  • fj.3  “ V fj.«  ■ V fj,5  ■ V f],6  ■ *j  <12) 


The  grid  pointe  in  physical  (x,y)  space  are  mapped  onto  a etraight  line  in 
s space.  Let  each  grid  point  be  connected  to  its  neighbor  in  s apace  by  a 
epring  with  spring  constant  between  the  j and  j+1  points  to  be 

defined  later. 
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Let  Anj  be  the  distance  between  the  j and  j+1  points  in  s apace. 

Then  Kj  An^  equals  a constant.  Therefore, 

An.  - Kx  An1/K.  (13) 

The  known  total  length  of  s can  be  expressed  as 

JN  JN 

sjh+i  * 2 ■ ^iki  2 i/Ki  “ AniKi  suk  ti4) 

j»l  j-1 

where  JN  ■ total  number  of  intervals.  Coneequently , 

Ani  • sjn*i/(ki  sum> 

Anj  » S^^tKj  SUM)  (15) 

The  constants  can  be  defined  explicitly  or  implicitly.  An  explicit 

definition  of  Kj  assigns  a value  to  the  spring  constant  independent  of  the 
position  of  the  epring  in  e space.  For  example, 

K,  ■ exp(-c^  • j/JN)  (16) 

provides  an  exponantially  increasing  spacing  of  grid  points  from  j * 1 to 
j » JN  where  the  constant  c^  is  used  to  control  stretching.  Thie  definition 
provides  no  adaptive  capability  but  dose  provide  a very  quick  way  of  con- 
centrating points  near  boundaries. 

An  implicit  definition  of  assigns  a value  to  the  spring  constant 

which  ie  a function  of  the  position  of  the  spring  in  s space.  For  an 
isiplicit  case,  the  algorithm  is  implemented  ae  followe: 

(1)  Start  at  the  first  column  (or  row)  of  data  and  compute  and  etora  all 
the  values  of  e.  and  f, 

(2)  Initialise  n and  n.  equal  to  sH  for  all  j,  ITER  « 0. 

J J 

(3)  Sat  ITER  - ITER  > 1. 

(4)  Using  the  known  stored  values  of  independent  variable  Sj  and  dependent 

variables  f,  , compute  the  new  values  of  the  dependent  variables  d.  at 
)«■  3,m 

ell  n.  with  a univariate  interpolation  routine  where  d . - x. , d _ * y , 

J j*ijj**j 

etc. 

(5)  Compute  K * K4(d,  ). 

J 3 J *■ 


028 


6)  Compute  the  new  coordinatee  nj  from  Equation  (15)  > set  “ e_j+j 

to  keep  end  point  valuee  unchanged. 

(7)  Compute  an  error  norm  to  determine  convergence 


L 


2 


("4 


(17) 


If  Lj  < 0.0001,  proceed  to  etep  (9). 

If  ITER  > 100,  stop  the  procedure. 

(S)  Set  n^  » n^  for  all  j and  go  to  step  (3). 

(9)  Replace  x. , y.,  p.,  etc.,  with  d.  , d.  ,,  d , etc.,  go  to 
J J J J*1  J « * 

next  column  of  data  and  repeat  all  etepe. 

Typical  definitione  of  employ  gradients  of  velocity,  internal  energy, 

or  Mach  number.  For  example. 


Kj  ” 1 + c3letnj+i>  " e(njH/(nj+1  ~ "j)  U8) 

where  c^  can  be  ueed  to  control  the  concentration  of  mesh  points  in  region.* 
of  high  gradient.  Examples  of  grid  dietribution  using  these  gradient  adaptive 
mechanisms  will  be  presented  in  the  RESULTS  section. 

In  caees  where  a high  gradient  develops  rapidly,  ae  in  the  vicinity  of  an 
expaneion  corner  of  a planetary  probe,  the  adaption  algorithm  can  have  trouble 
converging.  Thie  problem  is  overcome  as  follows. 

(1)  After  step  (4)  in  the  adaption  algorithm,  the  overall  change  in  the 
new  grid  distribution  n^  is  damped  by  writing 

n^  * c4n.  ♦ (1  - c4)n^  (19) 

where  0 < c4  £ 1 (e4  ie  typically  of  order  1/2) 

(2)  After  etep  (5)  in  the  adaption  algorithm,  the  values  of  the  spring 
constants  are  filtered  by  writing 

Kj"+1  “ <Kj-’  * 2Kj  * Kj»l)r'/4'  > “ 2 , JN-1  (20) 

Often  the  convergence  problems  can  be  overcome  by  keeping  the  reconing  interval 
small  or  by  avoiding  large  changee  of  c^  in  Equation  (18)  between  rcconings. 

The  major  advantage  of  the  rezoning  procedure  deecribed  herein  is  that  it 
provides  a very  quick  way  of  providing  an  adaptive  capability  to  the 
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, lation  The  option  to  control  zoning  frequency  allows  the 
finite-volume  formula  • featUres  of  the  flow  <!•*•>  3hock< 

9rid  to  dynamically  adapt  to  developinq  ^ ^ ^ ^ #pproached.  The 

ehear  layers)  and  then  to  he  • ^ ^ ,eparate  program  to  etudy  the  effecte 

re toning  procedure  can  different  definitions 

of  monitoring  various  gradiente  in  the  field 

of  K.. 

RESULTS  and  DISCUSSION  alqorithn  described  in  the  previous  sectione 

Tha  finite-volume,  adaptive  over  sphar„,  cylinders  and  two 

- - - *>  — TnrsSOnl Tst  cc^raheneive  teste  to  date  for 
plane trry  probe  configurat  on  ^rcshell  (Fig.  6)  . » 

teeting  the  adaption  algorithm  have  been  «*•  ^ ^ ^ ^ ^ 

configuration  which  lil^deal  exclusively  with  those  reeulte. 

Martian  atmoephere.  Th  - 1 285,  Re.  - 5000 ’ h#v*  b*,n  ch°Sen 

The  freestream  conditione,  “ 2'  Y ' ‘ “4.  Subeequent  calculatione 

results  of  Reference  i*- 

.nit  coep.rl.cn  -l«  th  a,TO1.trot.  the  cepehllltio. 

„ hl,».r  WnoU*  ■— »*  **“ 

*■  "TZlTZ* . - — — - ->“*  * "> 

’ T /n H 15  9.  point,  non.  th.  holy  ».  *"“>  -* 

116,1  "1Ul  C,/  Z hoOy/.eke  centerline  end  the  bo.  ebook. 

“ -yU  — 

The  preseure  dietributi  auDroximately  15  minute*  of 

....  was  run  to  convergence  in  approx i mat.  r 

Figure  7.  Thie  lQ-5  #ec/it.retion/grid  point.  Recent 

computing  time  with  a t ming  - stored  have  de- 

chengee  in  the  de  .««=««  *»  “ , ‘ ' 

creeeed  tdn,  to  .1..  * l"'5  „ ... 

in  el.pt IV,  gri.  deletion  -1th  W 

to  thi.  ~~  Ptobl-  In  the  ,n—  ^ " ^ 

— - - — - r::r; 

gince  th.  ^ „ run  thi.  =.e.  -1th  « **•“  “““  '*>' 

' - — - - v ::.:;.:r:c~.t^;: 
reVnUdUled^  , 

r «—  - — — - - 


830 


Reference  14  for  shock  shape  and  prsssurs  distribution  down  th  v 
show  gensrally  good  agreement  i *■  own  ths  wake  centerline 

y good  agresment  in  Figurss  9 and  7,  respectivsly 

Solutions  for  the  higher  Reynolds  number  cases  (5  • i0<  5'„  .-«. 

generated  ssguentially  using  the  previous  results  a.  initial  conditil 

C.U,  in  .rror  tbgy  u.  „„„  ro„in>,,  «“  *»  »»‘- 

g'Mimt.  .-.oci.t.a  with  „ch  flovJ  ' ° r”°1'''  °»  U"H 

s«*  int.r,.ti„,  of  th.„  hl5h 

In  riw„  10  through  17.  n.  co„»„^,  ,tld  5how.  , 

Point,  in  thn  ,h.,r  Uy„  ^ ^ MgJ  « 

Pnttern  nround  thn  nnpnnslon  v.hlcl.  md  ■" 

contours  p„,lnt.d  „ *”*  **  ' 

i.  datlned  by  th.  phy.ic.l  loc.tion  of  nxt.nt  of  U,„.  flours. 

r i”i"t  ai~- — - » *na 

r - - — — - Slvsn  in  Hour.  „ J‘ 

™::r  ^ — — — * <». «”  ■ 

“*  regions  i„  th.  , 

non  corner,  tb.  c.pturnd  r.co^r.s.ion  ,hoc„  in  th.  lu  h TT 

gredisnt  .hodc-liic,  region  on  tb.  ,mtty  u . ' * 

recirculeting  ™iociti«,  tepidly  cb„„  “*  ’**“ 

i P y cn*"9e  from  supersonic  to  subsonic  Th<„ 

phenomenon  was  el.o  observed  in  the  results  of  Reference  14 

numerical  epproech.  The  high  gredient  free  shear  1 ’*ln9  * dlfferent 

ona-fourth  maximum  body  radius  beyond  th.  corn.r  after^hTh^  ' 

«Pidly  decrease  and  the  flow  turn,  toward  the  axis.  ^1  ^ 

3“St  bel°W  sh3«  l»yer  of  low  density  , 2 < p/p  < . “**  r^10n 

* - number  distribution  along  „ „ ^ 

r r;:  rrr.  - — - 

ZTj.zziizr,z'“  “h  r‘,i““ in  zzzi, 

adaption  process.  ‘Srjr*te  Ver5etUU'  “d  «n.itivity  of  the 
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solution  at  every  point.  The  formulation  is  second  order  accurate  when  constant 
time  increment  is  used  to  advance  the  solution.  However,  at  large  Reynolds 
numbers  (small  cell  volumes)  the  constant  fit  advancement  of  the  solution, 
using  one  value  of  At  that  is  stable  for  all  points,  is  impractical  because 
many  more  iterations  are  required  to  converge  the  solution.  Converged  solutions 
for  complete  blunt  body  flows  using  approximately  3000  mesh  points  can  be 
obtained  within  30  minutes  using  the  variable  fit  option  and  adaptive  grid 
rezoning  on  the  CYBER  203  computer. 

In  all  of  the  problems  considered  herein  adaption  along  only  one  coordinate 
line  is  quite  sufficient  for  the  purpose  of  moving  mesh  points  to  high 
gradient  regions  in  the  flow.  Restricting  grid  motion  along  one  coordinate 
direction  permits  application  of  a highly  efficient  implicit  grid  adaption 
procedure.  The  routine  can  be  implemented  at  any  desired  frequency  thus 
permitting  dynamic  adaption  or  cost  savings  by  adapting  only  after  large  incre- 
ments of  iteration  count.  Grid  point  concentration  in  the  boundary  layer  and 
free  shear  layer  have  demonstrated  the  success  and  versatility  of  the  adaption 
algorithm. 
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pig.  1 Schematic  ehowing  relation 
of  indexing  system,  phyeicel  (x,y) 
coordinates,  and  computational 
(C,h)  coordinatee. 


4 Parepactiva  view  of 
asure  on  (x,y>  plana  for 
ique  shock  calculation  on 
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Fig.  2 Honintarective  example  of 
grid  adapting  to  prescribed 
praaaura  diatribution  with  two 
degraaa  of  adaptive  freedom. 
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Tig.  3 Interactive  example  of 
grid  adapting  to  coeluted  praaaura 
diatribution  acroaa  oblique  ahock 
with  two  degree#  of  adaptive 
freedom. 
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Fig.  5 Scheawtic  ehowing  apring 
■vatem  for  one  degree  of  adaptive 
freedom  along  n coordinate  line. 
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Fig.  6 Grid  over  Viking  Aeroshell 
obtained  with  explicit  definition 
of  K j , Cj/JN-.IS. 


Fig.  8 Grid  over  Viking  Aeroshell 
obtained  with  implicit  definition 
of  K , , Re  -5000. 
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Fig.  7 Pressure  distribution 
eround  body/weke  centerline  of 
Viking  Aeroshell.  M^-2 , Re— -5000, 
Y-1.285. 
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Fig.  9 Bow  shock  over  Viking 

Aeroshell  for  M -2,  Re-5000 . 
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Fig.  12  Mach  number  contour* 

around  expaneion  comer  of  Viking 

Aeroshell,  M -2 , Re  -106. 
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Fig.  10  Grid  over  Viking  Aeroehell 
obtained  with  implicit  definition 
ot  K , Re  -106. 
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Fig.  11  Streamline*  around 
expaneion  comer  of  Viking 
Aeroehell,  M *2  , Re  *>106. 
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Fig.  13  Velocity  and  grid  *pecing 
dietribution  along  ri  coordinate 
line  ahead  of  expansion  comer, 
Re^-lo4  and  10®. 
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Fig.  14  contour  plot  of  internal 
energy  over  Viking  A.  roshell  for 
M -2,  Re  “106. 
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Fig.  15  Mach  number  and  grid 
specing  distribution  through 
sepereted  shear  layer  behind 
expansion  comer  along  f)  coordinate 
line. 


